Spin moment formation and reduced orbital polarization in LaNiOs/LaAlOs 

superlattice: LDA+U study 



(N 
O 



I 

o 
o 



>: 
^, 

in; 

o: 



Myung Joon Han and Michel van Veenendaal 
Department of Physics, Northern Illinois University, De Kalb, Illmots 60115, USA and 
Advanced Photon Source, Argonne National Laboratory, 
9700 South Cass Avenue, Argonne, Illinois 60439, USA 
(Dated: January 26, 2012) 

Density functional band calculations have been performed to study LaNiOs/LaAlOs superlattices. 
Motivated by recent experiments reporting the magnetic and metal-insulator phase transition as a 
function of LaNiOa layer thickness, we examined the electronic structure, magnetic properties, and 
orbital occupation depending on the number of LaNiOa layers. Calculations show that the magnetic 
phase is stabler than the nonmagnetic for finite and positive U values. The orbital polarization is 
significantly reduced by U even in the magnetic regions. The implications of the results are discussed 
in comparison to recent experimental and theoretical studies within the limitations of the LDA+f/ 
method. 

PACS numbers: 73.20.-r, 75.70.-i, 71.15.Mb 



I. INTRODUCTION 

Understanding transition metal oxides is of perpetual 
interest and importance in condensed matter physics and 
material science Recent advances in layer- by- layer 
growth techniques of heterostructures of transition metal 
compounds have created particular interest due to their 
great scientific and technological potential Exotic 
material phenomena that are clearly distinctive from the 
'normal' phases include interface superconductivity Q, 
magnetism I, i, charge [1, Q and orbital reconstruc- 
tion 0. 

One of the most intriguing classes of materials may 
be the nickelate superlattices A series of re- 

cent theoretical studies have created considerable inter- 
est, suggesting the heterostructuring-induced orbital po- 
larization and the possible high-Tc superconductivity in 
LaNiOa (LA0)/LaA103 (LAO) superlattice In 
this picture, the two degenerate Ni eg orbital states 
are split by a combination of translational symmetry 
breaking and on-site Coulomb interaction, leading to a 
cuprate-like band structure. Although this kind of the- 
oretical picture has been challenged by more recent dy- 
namical mean-field theory (DMFT) calculations based on 
a charge transfer model including oxygen states explic- 
itly [l5|, several experimental papers have found other 
interesting phenomena in this system. Boris and co- 
workers [lO| reported a metal to insulator transition as a 
function of LNO layer thickness: Even though the bulk 
LNO is a paramagnetic (PM) metal, the heterostructure 
(LNO)„i/(LAO)„ becomes insulating and magnetic if m 
is small, m < 2, while it remains PM and metallic when 
TO > 4 [10| . The insulating behavior was also observed by 
Freeland et al. jl^. However the detailed magnetic and 
electronic structure changes, as well as the other impor- 
tant physical quantities such as orbital polarization have 
not yet been clearly understood as a function of layer 
thickness, m. 

In this study, we performed a detailed first-principles 



analysis for the electronic structure, magnetism, and or- 
bital polarization of LNO/LAO superlattice using the 
LDA-)-L/ method [l8l - l20| . Since previous calculations 
were performed at the LDA ^ or DMFT level [H dl] 
and assumed the bulk-like PM phase, LDA-|-t/ calcula- 
tion can provide meaningful information especially re- 
garding the magnetism in this system. Total energy cal- 
culations showed that the LDA-|-[/ stabilizes ferromag- 
netic (FM) spin order along in-plane and out-of-plane 
direction. This result may indicate the existence of an- 
other ground state configuration in between to ^ 2 and 
TO ~ 4 superlattices. The orbital polarization is signif- 
icantly reduced by U which is in contrast to the sim- 
plified Hubbard- type model prediction but consis- 
tent with the extended charge-transfer model DMFT 
calculation [l^. The calculated Ni-d valence based on 
LDA-t-t/ supports the recently suggested picture for the 
metal/insulator phase diagram based on the d- valency 
plj . and demonstrates the importance of the double 
counting issue. These results are discussed in compar- 
ison to the recent theoretical and experimental studies. 



II. COMPUTATIONAL DETAILS 

For the band-structure calculations, we employed 
TrouUier-Martins type norm-conserving pseudopotential 
[2^ with a partial core correction and linear combina- 
tion of the localized pseudo-atomic orbitals (LCPAO) [2^ 
as a basis set. In this pseudo-potential generation, the 
semi-core 3p electrons for transition metal atoms were 
included as valence electrons in order to take into ac- 
count the contribution of the semi-core states to the elec- 
tronic structure. We adopted the local density approxi- 
mation (LDA) for exchange-correlation energy functional 
as parametrized by Perdew and Zunger (23 | , and used en- 
ergy cutoff of 400 Ry and k-grid of 12 x 12 x 6 per unit 
superlattice volume. The LDA-I-C/ functional is adapted 
from the formalism of Ref. ^ and Ref. [H [H]. The 
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FIG. 1: (Color online) (a) Calculated total energy of 
(LNO)i/(LAO)i with the in- plane spin ordering of FM and 
AFM as a function of U. (b) Calculated total energy 
of (LNO)2/(LAO)i with the out-of-plane spin ordering of 
FM and AFM as a function of U where the in-plane spins 
are set to be FM. (c) Calculated Ni magnetic moment for 
(LNO)i/(LAO)i and (LNO)2/(LAO)i as a function of U (FM 
order considered), (d) Average number of Ni-Cg electrons in 
(LNO)i/(LAO)i and (LN0)2/(LA0)i as a function of U. 

geometry relaxation has been performed with the force 
criterion of 10"'^ Hartree/Bohr. During the relaxation 
process, the in-plane lattice constant is fixed considering 
the substrate effect in the experimental situation. The 
orbital polarization P, defined as 

can be calculated by integrating the projected density- 
of-states (DOS) up to Fermi level. All the calculations 
were performed using the density functional theory code 
OpenMX [281]. 

III. RESULT AND DISCUSSION 

One natural and evidently important question raised 
from experimental studies on the magnetic moment for- 
mation in this heterostructure [13] is about the ground 
state spin structure and its moment size. Fig. [1] sum- 
marizes our results. Fig. [TJa) shows the calculated total 
energy of (LNO)i/(LAO)i superlattice as a function of 
U. While at U=0 eV both FM and antiferromagnetic 
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FIG. 2: (Color online) The suggested schematic phase dia- 
gram (based on the previous experiment and our calculation) 
of the possible ground state configuration of (LNO)m/ (LAO)„ 
as a function of m (see the text). 

(AFM) spins eventually converge to a PM solution, the 
magnetic ground states are stabilized at the finite U. As 
clearly seen in Fig.[lja), FM spin ordering is energetically 
favored within the LNO plane. The energy difference be- 
tween FM and AFM state is 3, 246, and 525 meV per 
(LNO)i/(LAO)i for f/=3, 5, and 7 eV, respectively. FM 
spin order is also favored along the out-of-plane direction 
as shown in Fig. Wijo) in which we present total energies 
of (LN0)2/(LA0)i; the out-of-plane spin orderings are 
set to be FM or AFM while the in-plane order is FM. The 
calculated FM-AFM energy difference is 60, 126, and 201 
meV per (LN0)2/(LA0)i for U=3, 5, 7 eV, respectively. 
The energy differences between FM and AFM spin struc- 
ture becomes larger as U increases in both cases of in- 
plane and out-of-plane ordering. That is, the on-site cor- 
relations stabilize FM spin ordering, and as U decreases, 
the two magnetic solutions becomes more close in their 
energies, and eventually converges to the PM phase at 
U=0. The magnetic moment is also dependent on U. 
Fig. [He) shows the calculated magnetic moment for FM 
case as a function of U. The moment increases as U in- 
creases as in the other typical correlated transition metal 
oxide materials [13, [lO, U^- It is noted that, as soon as 
U is turned on, the Ni moment becomes ^ 1 fiB {U=3 
eV), and further increases to be 1.2 and 1.4 fiB at U=5 
and 7 eV, respectively. 

Our calculation results may seem to suggest that the 
magnetic moment formed in the thin-LNO superlattice 
[lo| is ordered ferromagnetically. Since the spin polarized 
oxide heterostructure could be useful for the device ap- 
plications, there has been active research for finding the 
structure that produces FM spin order [l^-lsH . Therefore 
our result of a FM ground state in LNO /LAO may have 
a positive implication for such an application. However 
it should be noted that the muon spin rotation (/iSR) 
experiment by Boris et al. [lo| is not well interpreted in 
the long range FM ordering picture even though the /xSR 
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FIG. 3: (Color online) Ni-eg DOS of (LNO),„/(LAO)i super- 
lattice geometry with PM spin (C/=0): (a) m = 1, (b) m = 2, 
(c) m — 3, and (d) m = 4. Solid (red) and dotted (blue) 
lines correspond to d^2_y2 and d^^2_j.2 states, respectively, 
and Fermi level is set to be 0. 

is basically a local probe and that the origin of nietal- 
insulator phase transition in the nickelate series are not 
clearly understood yet. Especially regarding the charge 
disproportion or ordering in nickelates, the conventional 
LDA+U has a clear limitation to describe such phenom- 
ena [12, H^- Moreover Fig. [IJa) and (b) indicate the 
possibility of AFM ground state in the negative U region 
which is more or less related to the reported charge dis- 
orders in the nickelate systems [13, [13|- Therefore one 
needs to be careful in the interpretation of our LDA+U 
results on the FM spin ground state as an indication of 
long range ordered ground state as one may see in the 
actinide systems for example [sij, [ssj . 

Our results have another interesting implication re- 
garding the phase diagram. Since the thin-LNO superlat- 
tice, (LN0)m<2/(LA0)„, can have either FM insulating 
(FM-I) ground state as predicted by LDA+U calcula- 
tions, or, AFM insulating (AFM-I) one which is more 
consistent with the /iSR experiment, the system would 
have FM metalhc (FM-M) or AFM metalhc (AFM-M) 
region in between the thin-LNO limit (to < 2), and the 
bulk-like thick-LNO limit (m > 4; PM and metallic, PM- 
M) . One may also expect the PM insulating (PM-I) 
phase stabilized in the same region of m; the intermedi- 
ate regime, 2 < m < 4 (see. Fig. [2]). As LDA-I-J7 method 
is unable to describe correlated PM solutions properly 
[H, [13 , the further pursue along this line is beyond the 
scope of our study. For that purpose, one may resort the 
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FIG. 4: (Color online) Ni-e^ DOS of (LNO)^/(LAO)i su- 
perlattice geometry with FM spin ((7=5 eV): (a) m — 1, (b) 
m = 2, (c) m = 3, and (d) m = 4. Up/down panels repre- 
sent the up/down spin states, respectively. Solid (red) and 
dotted (blue) lines correspond to d^2_y2 and d^^2_^2 states, 
respectively, and the Fermi level is set to be 0. 



dynamical mean field theor y (D MFT) calculations with 
charge self-consistency [H,113], and compare the total 
energy for PM, FM and AFM configurations as we did 
in this study within LDA+U scheme. 

To understand metal-insulator phase transition as a 
function of LNO thickness [l^l, we examined Ni Cg DOS 
depending on the thickness, m. The d-band width change 
as a function of to may be an useful information as U/W 
plays an important role in the metal-insulator phase tran- 
sition of rare-earth nickelate Fig. [3] shows DOS of 
the PM case with J7 = 0. Even if the bulk LNO locates 
at the vicinity of metal-insulator phase boundary and 
exhibits the correlated electron behaviors, the electronic 
structure of PM LNO has been reasonably well described 
within LDA (or GGA) as shown in the previous studies 
[13, [H, [HI- Fig. |31[a)-(d) presents the evolution of eg 
DOS as a function of m. It is noted that the d^^2_j.2 
band width notably changes, whereas the (1^2 _y2 width 
remains almost same across to=2-4. The to = 2 case 
(Fig. E^b)), which was reported to be magnetic and in- 
sulating [10], the right edge of the c?3^2_^2 state is fairly 
similar to that of dx2-y2 whereas the left-edge is reduced. 
It is not certain however that such a relatively small dif- 
ference in terms of the effective band width is responsible 
for the metal-insulator phase transition between m = 2 
and 4 as reported in experiment [lo| . 
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FIG. 5: (Color online) The averaged orbital polarization as a 
function of LNO layer thickness. Dotted (blue, double cross), 
double-dotted (black, circle), and solid (red, square) lines rep- 
resent U=0 (PM), (7=3 (FM), and U=5 (FM), respectively. 



A complementary picture can be provided by LDA+U 
calculations for the magnetic phase. Fig. |4l[a)-(d) shows 
the evolution of Cg DOS as a function of m for the FM 
(LNO)i/ (LAO)i. Once again the notable change is found 
in d^z^^r^ band; for m = 1, the up-spin d^z^^^.^ state 
(upper panel) forms a fairly localized DOS around Fermi 
level. This state becomes more and more delocalized as m 
increases as seen in Fig. Ill^b)-(d). Once again, however, 
the amount of the effective band width change depending 
on m does not seem to be enough to make metal-insulator 
phase transition across m=2-4. Other origins than the 
simple U/W change may be more relevant to the transi- 
tion [2ll,l33|. 

An interesting point observed in Fig. 2] is that even in 
LDA+U calculation, the systems do not become a perfect 
insulator, but have finite number of states around Fermi 
level [i^l- It may partly be attributed to the strong co- 
valency between Ni-d and 0-p states. As 0-p is much 
less affected by U, the small amount of DOS is not re- 
moved perfectly. Another important factor is the double- 
counting energy correction for which several functional 
forms have been suggested |T8U23, [2^, [2^ , but there is no 
well defined solution still yet [2l|: As the double-counting 
energy is typically represented by ^UNd{Nd — 1), the 
LDA+U charge self-consistency adjusts Nd (or, the ef- 
fective d-level energy and therefore the charge transfer 
energy, ep — e^), depending on U. According to the re- 
cent DMFT study, which tunes the double-counting term 
as a parameter^the system eventually becomes insulat- 
ing at large U [la, |4l| . Therefore the small states around 
the Fermi level obtained by LDA+U can be attributed to 
the double-counting error that is hardly handled within 
the current formalism of LDA+U. 

An interesting recent finding on (3-dimensional) nick- 
elates is that the metal-insulator transition occurs at the 
very narrow region of Nd and the same holds for the 
cuprates [2l|. From the single-site DMFT calculations 



with a double-counting energy as another tuning param- 
eter, Wang et al. [2l[ showed that there is a well defined 



iVc, that is, the critical value of d-valency of transition 
metal; the system is metallic if Nd > N^ and insulat- 
ing if Nd < Nc. For nickelates, Nc" ~ 1.3. That is, 
for the parameters (implicitly including double-counting 
correction), which result in Nd > Nc, the system re- 
mains metallic even for the very large U. This conclu- 
sion suggests Nd as the critical variable for understand- 
ing charge-transfer systems [21i] . Now it might be in- 
structive to analyze our LDA+U results within this new 
picture of Ref. [2l|. Even if LDA-l-17 is a lower level ap- 
proximation compared to DMFT, a big merit of LDA-I-J7 
is that it can be performed with the whole charge self- 
consistency and take magnetism into account while the 
previous DMFT calculations have dealt with the PM 
phase [T^, HI]. Fig. md) presents Nd for PM {U=0) 
and FM {U=3, 5, 7) calculations of (LNO)™=i/(LAO)i 
and (LN0)m=2/(LA0)i. While Nd decreases as U in- 
creases, the difference is not significant; AA^^ w 0.15. 
Therefore the current implementation of LDA+U and 
the charge self-consistency roughly follow the constant 
Nd line. Now we note that the value of N^^ > 2.2 is far 
from the A^c « 1-3 predicted by DMFT. Therefore it is 
consistent with the metallic ground state and supports 
the conclusion of Ref. [2ll 

Relative orbital occupation in the two Cg orbitals is 
a central quantity in understanding LNO/LAO systems 
[1 [H H (Hi- The orbital polarization can be non-zero 
due to the translation symmetry breaking induced by 
heterostructuring while the bulk polarization is 0. Ac- 
cording to the previous calculations, highly polarized 
P can possibly drives the system to be high-Tc super- 
conductor jl^, On the other hand, a more recent 
DMFT calculation based on the realistic model hamil- 
tonian including oxygen orbitals predicted that the po- 
larization is actually reduced by the on-site correlation 
[l5| . Since both of the previous DMFT calculations as- 
sumed PM phase and did not consider the full charge 
self-consistency, LDA+U calculation can give a compli- 
mentary information for the magnetic solution in spite 
of its limitation in describing correlation effect compared 
to DMFT. The result is presented in Fig. [51 It is noted 
that the inclusion of U reduces the polarization signif- 
icantly. There is a large separation between U=0 and 
U—3 results, while the differences between U=3 and 5 
results are small. Importantly the calculated polariza- 
tions for finite U are order-of-magnitude same compared 
with the recent DMFT results [15]. Therefore it supports 
the conclusion of the recent DMFT calculation [I^ . The 
small polarization at the finite U can also be seen in DOS 
presented in Fig. 2] where no big difference can be found 
between the two eg orbital occupations. We note that 
the reduced orbital polarization by U can be compati- 
ble with the ferromagnetic spin order. For cuprates, for 
example, the fully polarized — j/^ orbital is directly re- 
lated to the in-plane AFM spin order. In the nickelates, 
we have one more orbital degree of freedom available and 
the spin can align ferromagnetically. And actually it is 
found that the enhanced ferromagnetic trend by increas- 
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ing U corresponds to the reduced orbital polarization as 
shown in Fig. [5j 

IV. SUMMARY 

Using the band structure calculations based on 
LDA+C/, we examined the magnetic moment, electronic 
structure, d-valence, and orbital polarization as a func- 
tion of U and m (LNO thickness). The calculated results 
clearly showed the formation of magnetic moment at the 
finite U region being consistent with a recent /xSR, but 
the long range ordering pattern is not so clear consid- 
ering the experiment. While d^^^-r^ band width is re- 
duced as m approaches to 1, it may not be enough to 



be responsible for the metal-insulator transition. The 
calculated orbital polarization is significantly reduced by 
U, strongly supporting the conclusion of a recent DMFT 
calculations and indicating the absence of high tempera- 
ture superconductivity in this system. 
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